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ABSTRACT 

We estimate the optical depth to self-lensing by stars in the Galactic bulge 
using the Hubble Space Telescope star counts of Holtzman et al. and Zoccali 
et al. as extrapolated by Gould into the brown-dwarf and remnant regimes 
and deprojected along the line of sight using the model of Dwek et al. We 
find a self-lensing optical depth r(bulge — bulge) = 0.98 x 1CT 6 . When com- 
bined with the lensing of bulge stars by foreground stars in the disk, this yields 
r(bulge — total) = 1.63 x 10~ 6 , in reasonable agreement with the estimates of 
r = 2. 13 ±0.40 x 10~ 6 and r = 1.08 ±0.30 x 10~ 6 based on observations of clump 
giants by the MACHO and EROS collaborations. 

Subject headings: gravitational lensing - stars: luminosity function, mass func- 
tion Galaxy: bulge - dark matter 



1. Introduction 

The optical depth to microlensing toward the Galactic bulge has been controversial since 
before it was first "officially" measured. When Paczyhski (1991) and Griest et al. (1991) 
first proposed bulge microlensing observations, they estimated the optical depth to be about 
r~5x 10~ 7 assuming all events were due to known disk stars. However, as soon as the first 
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six bulge events were reported by OGLE (Udalski et al. 1993, 1994), the apparently high 
event rate prompted Kiraga & Paczyhski (1994) to evaluate the contribution to the optical 
depth of bulge stars in addition to those of the disk. They estimated r ~ 8.5 x 10~ 7 and 
concluded that the value could be as much as twice as high if the bulge were elongated along 
the line of sight. 

Nevertheless, the first measurements of the optical depth, r = 3.3 ± 1.2 x 10~ 6 by 
OGLE (Udalski et al. 1994), and r = 3.9+{| x 1(T 6 by MACHO (Alcock et al. 1997) 1 were 
substantially higher than even the highest predictions. Indeed, Gould (1994) and Kuijken 
(1997) argued that they were so high as to eliminate the need for any dark halo inside the solar 
circle, provided the Galaxy was assumed axisymmetric. This argument was subsequently 
generalized to non-axisymmetric mass distributions by Binney, Bissantz & Gerhard (2000). 

These early results motivated substantial additional work. On the one hand MACHO 
made several new estimates of the bulge optical depth, one based on a difference-image 
analysis of a large number of primarily faint sources, yielding r = 3.2 ± 0.5 x 10~ 6 (Alcock 
et al. 2000), and another based on clump giants, yielding r = 2.0 ± 0.4 x 10~ 6 (Popowski et 
al. 2001). On the other hand, beginning with Zhao, Spergel & Rich (1995), many workers 
attempted to develop highly non-axisymmetric bulge models that could account for the high 
observed optical depth (Metcalf 1995; Zhao, Rich & Spergel 1996; Zhao & Mao 1996; Bissantz 
et al. 1997; Gyuk 1999; Nair & Miralda-Escude 1999; Sevenstar & Kalnajs 2001; Belokurov 
& Evans 2002). 

Most recently, EROS has reported a much lower optical depth, r = 0.94 ± 0.30 x 10~ 6 
(Afonso et al. 2003) based on a study specifically designed to monitor clump giant sources 
as advocated by Gould (1995). In their Table 2, Afonso et al. (2003) also put all optical- 
depth measurements on a common basis by adjusting them for their offset from Baade's 
Window (BW). Hence, for example, their own measurement was adjusted upward to r = 
1.08 ± 0.30 x 10~ 6 because their mean field position was less densely populated than BW. 

Given the strongly divergent observational results as well as the theoretical difficulties 
these pose, it is useful to have as many model-independent constraints as possible. One 
such constraint comes from star counts. At the time that Paczyhski (1991), Griest et al. 
(1991), and Kiraga & Paczyhski (1994) made their estimates, the stellar content of the disk 
was poorly measured and the stellar content of the bulge was almost completely uncon- 
strained. Star count from Hubble Space Telescope (HST) observations have now completely 
transformed that situation (Holtzman et al. 1998; Zoccali et al. 2000; Zheng et al. 2002). 



1 Submission of the manuscript was delayed 4 months by a secretarial error (C. Alcock 1996, private 
communication) but the results were widely circulated in the microlensing community (astro-ph/9512146). 
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The stellar mass function is well measured in the disk down to the hydrogen-burning limit 
and in the bulge down to M ~ 0.15 M . There are, of course, uncertainties in extrapolating 
these mass functions to lower, sub-stellar masses, and to the remnants of the now deceased 
stars at higher masses. Nevertheless, the fact that large portions of these mass functions 
have been directly measured allows one to construct useful constraints. 

We formulate these constraints within the context of a specific model of the Galactic 
bulge and disk. The model of the disk is relatively secure, and uncertainties in it play an 
overall very small role in controversies about the total optical depth toward the bulge. That 
is why most of the effort to explain the high optical depth has centered on models of the 
bulge. The benchmark model that we adopt is therefore far from unique. We show, however, 
that it is possible to factor the stellar constraint into two terms, one representing the bulge 
model convolved with the observational strategy, and the other representing the star counts. 
In this way, our result can easily be applied to any model of the Galactic bulge. 



2. Bulge and Disk Models 
2.1. Disk 

We model the local vertical disk density profile in accord with the model of Zheng et al. 
(2002). To extend this model to the whole Galactic disk, we assume that the column density 
of the disk has a scale length H = 2.75 kpc, as measured by Zheng et al. (2002). We account 
for the gradual flaring of the disk by rescaling the scale heights in the Zheng et al. (2002) 
formula in proportion to the scale height derived by Kent (1992). We normalize the local 
stellar column density to S = 36M pc~ 2 . This includes 28M pc" 2 in observable stars 
and white dwarfs (Zheng et al. 2002; Gould, Bahcall, & Flynn 1997) and another 8 M Q pc -2 , 
which is a rough estimate of the column density of brown dwarfs (BDs). The disk density 
profile in cylindrical coordinates is then, 

P{R , Z) = ^ exp (_^)[ (1 _^ ech ^ + i3exp (_M)] (1) 



where p = 0.0493 M pc" 3 , (3 = 0.565, h x = 270 pc, h 2 = 440 pc, H = 2.75 kpc, R = 8 kpc, 
and 

'' W = maX {9ol^ + 114 - 670 }' (2) 
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2.2. Bulge 

For the bulge model, we follow Han & Gould (1995) and scale the bulge mass density 
to the deprojected infrared light density profile of Dwek et al. (1995). Specifically, we use 
model G2 (with -R max = 5kpc) from their Table 2. However, whereas Han & Gould (1995) 
fixed the normalization of this model using the tensor virial theorem (which is sensitive to 
all the mass), we will normalize it by the observed stars (and inferred stellar remnants and 
substellar objects). 

2.3. Mass Function and Mass-Luminosity Relation 

For the bulge, we adopt the (unnormalized) mass function of Gould (2000). This assumes 
that bulge stars formed initially according to a power law, dN/dM = k(M/M hrk ) a , where 
M hrk = 0.7 M Q , a = -2.0 for M > M brk , and a = -1.3 for M < M brk . These slopes are 
consistent with the observations of Zoccali et al. (2000), but the profile is extended below 
their lower limit of M ~ 0.15 M to a BD cutoff of M = 0.03 M Q . The stars with initial 
masses 1M < M < 8 M Q are assumed to have become white dwarfs (WDs), those with 
masses 8 M & < M < 40 M are assumed to have become neutron stars (NSs) with masses 
M = 1.35 Mq, and those with masses M > 40 M© are assumed to have become black holes 
(BHs) with mass M = 5M . The total mass is then dominated by main-sequence (MS) 
stars, with mass fractions 

BD : MS : WD : NS : BH = 7 : 62 : 22 : 6 : 3. (3) 

This is important because, as we will see in § 3, the observational constraint comes from 
luminous (i.e. MS) stars. 

We adopt the mass- My relation of Cox (1999) 

3. Normalization From Star Counts 

We populate the bulge with stars (and BDs and remnants) according to the Gould 
(2000) mass function described in § 2.3 and adjust the overall normalization until we match 
the Holtzman et al. (1998) V^-band star counts in BW. See Figure 1. To predict these counts, 
we assign each star a luminosity using the Cox (1999) mass- My relation, while treating all 
BDs, WDs, NSs, and BHs as dark. We then convert to apparent magnitudes using the each 
star's individual distance. Finally, each star is then reddened by assuming that the total 
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extinction along the line of sight is Ay = 1.28 (Holtzman et al. 1998), and that the dust has 
a scale height of 120 pc. 

We incorporate disk as well as bulge stars in predicting the Holtzman et al. (1998) star 
counts. Since the disk stellar profile is regarded as well measured (by Zheng et al. 2002), we 
do not adjust the normalization of the disk profile as we do the bulge profile, but rather leave 
it fixed in the form given in § 2.1. However, we use the same mass function and mass- My 
relation for the disk as the bulge. In principle, one should make an independent estimate 
of these functions. However, since the disk stars contribute only ~ 15% of the counts (see 
Fig. 1), the net corrections from more accurate functions would be only a few percent, which 
is small compared to other uncertainties in the problem. Hence, we ignore this distinction 
in the interest of simplicity. 

We fix the normalization by demanding agreement between the model predictions and 
the (mass-weighted) Holtzman et al. (1998) star counts over the range 22.5 < V < 26.5. At 
fainter magnitudes, the Holtzman et al. (1998) data become seriously incomplete, while at 
somewhat brighter magnitudes our simple mass-My relation fails to account for evolution 
off the MS and so slightly overpredicts the counts. At much brighter magnitudes, the fact 
that our model has no giants causes it to completely underestimate the counts. These latter 
two effects are each small and roughly cancel one another. See Figure 1. 

The good agreement over four magnitudes demonstrates that the mass function of Zoc- 
cali et al. (2000), which is based on infrared observations, is compatible with the optically- 
based Holtzman et al. (1998) mass function, an agreement already noted by Zoccali et al. 
(2000). 

The total column density of bulge stars (and associated BDs and remnants) toward BW 

is 

£* = 2086 M pc" 2 . (4) 

As shown by equation (3), 62% of this mass is in the form of luminous stars, while 55% is in 
MS stars with masses 0.15 M < M < 1 M , and so with sufficient luminosity to have been 
directly observed by Zoccali et al. (2000). 

4. Optical Depth Due to Stars 

The observed optical depth will always be an average over the individual optical depths 
to the stars being monitored. More distant stars have higher optical depth, but are also 
fainter and so less likely to be included in the sample (Kiraga & Paczyhski 1994). To 
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parameterize this effect, we write 

( v . 4ttG / °° dD s D^p(D s ) f D ° dD lP {Di)Di{D a - Di)/D a 
c 2 f Q d£> 8 £> s 1 p{D s ) 

For standard candles, which can be identified independent of distance, 7 = 0. We are most 
directly interested in comparing to optical depth measurements using clump giants, which 
are approximately standard candles. We therefore adopt 7 = as the default. For bulge 
sources and for, respectively, disk and bulge lenses, we find 

(r)o(disk) = 0.65 x 1(T 6 , (r) (bulge) = 0.98 x 10~ 6 . (6) 

In fact, (r) 7 (disk) does not significantly depend on 7, but (r) 7 (bulge) does. For example 
(r)i(bulge) = 0.86 x lO" 6 

The total optical depth due to stars that is predicted by this mass profile for observations 
toward BW, (r) (total) = 1.63 x 10 -6 , agrees reasonably well with the two measurements 
made using clump giants. When these are adjusted to a common mean direction at BW 
(Afonso et al. 2003), they yield r = 2.13 ± 0.40 x 10~ 6 (Popowski et al. 2001) and r = 
1.08 ± 0.30 x 10- 6 (Afonso et al. 2003). 



5. Discussion 

While the agreement between the observations and the model predictions is comforting, 
it is important to recognize that the model has some uncertainties. These are most easily 
discussed by writing the predicted optical depth as 



(r) 7 (bulge) = -JLJL. (7) 



This equation serves to define D T , which may be thought of roughly as the characteristic 
source-lens separation. This characteristic separation depends only on the mass profile and 
(through 7) on the observational strategy. It can be evaluated trivially for any mass model. 
For the Dwek et al. (1995) bulge model used here, 

A) = 782pc. (8) 



On the other hand, depends only on the surface density of stars. For luminous stars, 
this is practically an observed quantity (Holtzman et al. 1998; Zoccali et al. 2000), and the 
major uncertainty is how to extrapolate the observations into the BD and remnant regimes. 
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One might, for example, argue that a Salpeter mass function (a = 2.35) is more appropriate 
for the extrapolation to higher masses. In this case, the WD contribution would fall by 20% 
and the BH contribution would fall by a factor 5. Or one could argue that the mean mass 
of BHs is higher, or that the threshold stellar mass to leave BH remnants is lower than the 
40 M & adopted here, either of which would raise the BH contribution. Finally, the Zoccali et 
al. (2000) mass function does not take account of binary companions, which might plausibly 
raise by 10-20%. Any of these changes can be easily incorporated using equation (3), 
once new estimates are made. 

We can compare the approach adopted here to that of Han & Gould (1995), who nor- 
malized the Dwek et al. (1995) model using the tensor virial theorem. To do so, we first note 
that Han & Gould (1995) tilted their bulge ellipsoid 20° to the line of sight, as advocated in 
the Dwek et al. (1995) abstract, rather than the 13.°4 given in Table 1 of Dwek et al. (1995) 
and used here. Normalizing the Han & Gould (1995) 20° model according to the procedure 
defined in this paper, we find a bulge self-lensing optical depth of (r)o = 0.88 x 10~ 6 , only 
66% of the value obtained using the Han & Gould (1995) tensor virial theorem normalization. 
Part of the difference is undoubtedly due to binary companions, which contribute to the inte- 
grated mass and enter the microlensing optical depth but which, because they are generally 
substantially less luminous than their primaries, leave almost no trace on the observed HST 
stellar luminosity functions. If the BD or remnant populations were substantially larger than 
modeled by Gould (2000), this would also add to both the integrated mass density and the 
optical depth. Part of the difference could also be due to non-baryonic dark matter, which 
would add to the integrated mass (to which the tensor virial theorem is sensitive) but would 
not contribute to microlensing. Finally, of course, the tensor virial theorem yields results 
that are only as accurate as the mass model and velocity dispersion tensor to which it is 
applied, and either of these could have errors. 

Work by CH was supported by the Astrophysical Research Center for the Structure and 
Evolution of the Cosmos (ARCSEC) of Korea Science and Engineering Foundation (KOSEF) 
through the Science Research Center (SRC) program. Work by AG was supported by grant 
AST 02-01266 from the NSF. 
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Fig. 1. — Star counts toward Baade's Window. Observations of Holtzman et al. (1998) 
(connected data points) are used to fix the normalization of the bulge mass-density profile of 
Dwek et al. (1995). The individual contributions of the disk (dashed curve) and bulge (solid 
curve) are added to predict the total disk+bulge counts (bold curve). The disk is held fixed, 
while the bulge is scaled until the mass-weighted disk+bulge star-count predictions agree 
with the observations over the range 22.5 < V < 26.5. 



